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ABSTRACT We have produced two-dimensional maps of the inten- 
sity ratio Qqq of 60 /mi infrared to 20 cm radio continuum emission, for 
a set of 25 nearby galaxies, mostly spirals. The ratio maps were obtained 
from infrared images made using Infrared Astronomical Satellite (IRAS) 
data with the Maximum Correlation Method, and radio images made us- 
ing VLA data. Before taking the ratio, the radio images were processed 
so as to have the same resolution properties as the infrared images; the 
final spatial resolution in all cases is approximately 1', corresponding to 
1-2 kpc for most galaxies. These images allow us to study the variations 
of the Qeo ratio with unprecedented spatial resolution, and thus represent 
a major improvement over earlier work. 

Our new high-resolution maps confirm the slow decrease of Q $ o with 
increasing radial distance from the nucleus, but show additional structure 
which is probably associated with separate sites of active star formation in 
the spiral arms. The maps show Qeo to be more closely related to infrared 
surface brightness than to the radial distance in the galaxy disk. We 
expect that the results will provide improved constraints on the evolution 
(diffusion, decay and escape) of cosmic-ray electrons in the magnetic field 
of the disks. 


INTRODUCTION 

The strong correlation, observed in spiral galaxies, between flux densities at far- 
infrared and centimeter wavelengths (Helou, Soifer, and Rowan- Robinson 1985) 
has prompted comparisons between the spatial structure in these two wavelength 
regimes. In a recent study, Bicay and Helou (1990) compared the spatial distri- 
butions, at 60 fin i and 20 cm wavelength, for 25 galaxies (24 late-type spirals and 
one irregular). They found that the radio images had the appearance of smeared 
versions of the infrared images, and that the smearing could be understood in 
terms of the diffusion of cosmic ray electrons. In their model, discussed in more 
detail by Helou and Bicay (1993), the correlation between infrared and radio 
emission is driven by massive stars heating up the dust, and causing cosmic-ray 
electrons to be accelerated during the supernova phase; the cosmic rays then give 
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rise to synchrotron radio emission which dominates at 20 cm. The difference in 
spatial distribution between infrared and radio arises because cosmic rays diffuse 
up to a few kpc before being lost to radiative decay or escape. 

Bicay and Helou (1990) used the raw detector data from the Infrared As- 
tronomy Satellite (IRAS), which are of rather coarse spatial resolution (1.5' X 
4.7'). Since the latter dimension was comparable to the sizes of the galaxies 
themselves, the study was basically one- dimensional in terms of spatial struc- 
ture. In this paper we describe an improved approach to the same question 
using infrared data of higher spatial resolution (~ T in both coordinates). 

DATA 

The galaxy sample used in the current investigation consisted of the same 25 
galaxies studied by Bicay and Helou (1990). Our input data for each galaxy 
consisted of: 

(1) Radio continuum image at 20 cm wavelength (frequency 1.49 GHz) made 
using the VLA, and described by Condon (1987). The CLEAN algorithm was 
used for deconvolution, with circular Gaussian restoring beams of FWHM 0.8', 
0.9', or 1.0'. 

(2) Infrared continuum image at 60 gm wavelength, made using data from 
the Infrared Astronomy Satellite (IRAS), and the imaging technique known as 
the Maximum Correlation Method (MCM), described by Aumann, Fowler, and 
Melnyk (1990). The number of iterations used was 20. 

Although the infrared and radio images had comparable spatial resolution 
(~ 1'), the infrared images were nonisoplanatic. The behavior of their resulting 
point spread functions (PSF) was determined by the IRAS scanning geometry, 
the spatial responses of the set of IRAS detectors, and the inherent properties of 
the MCM. Since the goal of the study was to compare radio and infrared images 
at the same spatial resolution, further processing was necessary. 


ANALYSIS PROCEDURE 

In order to produce pairs of images (radio and infrared) with the same spatial 
resolution, we re-CLEANed the radio images and restored them using an effective 
PSF corresponding to the infrared images. The procedure involved the following 
steps: 

(1) Deconvolve the radio PSF. This step was necessary since the original 
sets of CLEAN components involved in Condon’s (1987) maps were not avail- 
able. The deconvolution was performed by a reapplication of the CLEAN algo- 
rithm, based on the same circular Gaussian beam which Condon had used to 
restore each map. The image was CLEANed down to the residual noise level 
(0.1 m Jy/ beam), yielding a source model in the form of a set of delta- function 
components. 

(2) Generate synthetic “raw detector data” as would be observed if this 
source model were scanned using detectors whose spatial responses were the 
same as for IRAS, and using the same scanning geometry as for the original 
IRAS observations. 
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(3) Run these data through the same MCM imaging algorithm as was used 
to generate the infrared images. The pixel size (15"), field of view, and position 
of the center pixel, were the same as for the corresponding infrared image. 

The result in each case was a radio image with spatial resolution properties 
similar to that of the infrared image. 

Because of uncertainties in the absolute positioning of the IRAS data (re- 
sulting in errors comparable to a pixel width in the infrared images), it was 
necessary to perform a final registration step. The required position offset of 
the infrared image, (Aa,A<5), was estimated by minimizing the sum of squares 
of residuals, </>(Aa, A<5, q), defined by: 

<f){Aa,A6,q) = ]T][/ 20 cm(«, <5) - 9 _ U 60 (a - Aa, 6 - A6)] 2 (1) 

QyS 

where /20cm( a ,<5) and I 60 (a,6) represent the intensity distributions at wave- 
lengths 20 cm and 60 /im, respectively, as a function of right ascension, a, and 
declination, <5, and q is a scaling factor representing the weighted-mean ratio 
of infrared to radio intensity. The summation was performed over all positions 
which fell within a rectangular box, constructed so as to include the galaxy but 
exclude background sources. 

The final numerical step was to calculate the infrared-to-radio intensity 
ratio, Qgo = ho/hocm , as a function of position. This quantity was then plotted, 
in the form of a greyscale image, for all positions for which Qeo > 3a q, where 
uq is the standard deviation in Q. The quantity oq was determined by the noise 
levels in the infrared and radio images, which were dominated by the dynamic 
ranges of the processed radio maps, typically 70:1. 

RESULTS 

Figure 1 shows the results for three sample cases, namely NGC 5194, NGC 5236, 
and NGC 6946. 

We have compared all 25 ratio maps with the corresponding infrared inten- 
sity maps, and find that in nearly all cases, there is a local maximum in Qeo 
coincident with the galactic nucleus, and that the average o in an annulus 
centered on the nucleus decreases monotonically with radial distance, 7‘, from 
the nucleus, in a manner which is in qualitative agreement with the Bicay and 
Helou (1990) model. The principal exceptions to this behavior were IC 10 (an 
irregular galaxy) and NGC 3031 (a disk galaxy with a strong Seyfert nucleus), 
neither of which were expected to fit the model. 

In addition to the general radial behavior of Qeo described above, the im- 
proved spatial resolution reveals in some galaxies off-nucleus local maxima in 
Qqq , which coincide with local maxima in /^o- These new details are seen in 
the well-resolved face-on spirals, and are particularly prominent in NGC 5194, 
NGC 5236, and NGC 6946. They also appear in the edge-on spiral, NGC 4631. 
These local maxima correspond to locally prominent star formation regions in 
the spiral arms. The new features cannot be accounted for in the simple model 
by Helou and Bicay (1993), which describes the / 6 o distribution as an exponen- 
tial disk. They therefore partially destroy the predicted gradient of Q 60 as a 
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FIGURE 1 Maps of 60 pm infrared surface brightness, presented on a loga- 
rithmic scale of 2.5 decades (left hand plots) and the intensity ratio of 60 pm 
infrared to 20 cm radio emission, presented on a Unear scale (right hand plots). 
The field of view in aU cases is 16' X 16'. (a) NGC 5194 (Peak value of ratio map 
= 284); (b) NGC 5236 (Peak value of ratio map = 266); (c) NGC 6946 (Peak 
value of ratio map = 220). 
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FIGURE 2 Scatter plots of (^0 as a function of radial distance from nucleus, 
?*, and infrared surface brightness, 1$ o, for NGC 5194, NGC 5236, and NGC 
69 16. 



120 


MARSH AND HELOU 


function of r. This is clearly visible in Figure 2, which shows scatter plots of 
Q 60 w. r for NGC 5194, 5236 and 6946. 

The essential content, however, of the Bicay and Helou model is not in 
the geometry, but in the physics of diffusion, escape and decay of cosmic-ray 
electrons, and the effects on the appearance of the radio map. These effects 
should apply locally to each center of star formation activity, and cause the 
ratio Q$q to fall with increasing distance from each center. Since / 6 0 also falls 
with increasing distance from a center of activity, one would expect <3 so and / 6 0 
to be positively correlated. We test this hypothesis using scatter plots of Q 60 
v. r and Q $ 0 v. / 60 . The results for NGC 5194, NGC 5236, and NGC 6946 are 
shown in Figure 2. The three plots are clearly consistent with our expectations, 
i.e., Qi 5 o is more closely related to infrared surface brightness than to the radial 
distance in the galaxy disk. 

In a forthcoming paper we will present the complete results for all 25 galax- 
ies in the sample, and discuss the physical constraints that these results impose 
on the evolution of cosmic rays in the magnetic field of the disks. 
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